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Abstract 

Unlike a helicopter, an insect can, in theory, use both lift and drag to stay aloft. Here we 
show that a dragonfly uses mostly drag to hover by employing asymmetric up and down strokes. 
Computations of a family of strokes further show that using drag can be as efficient as using lift at 
the low Reynolds number regime appropriate for insects. Previously, asymmetric strokes employed 
by dragonflies were viewed as an exception. Here we suggest that these asymmetric strokes are 
building blocks of the commonly seen symmetric strokes, such as a figure-of-eight or a U-shape. 
Therefore insects which use those symmetric strokes can also use some drag to hover. In a sense, 
some insects are rowers or swimmers in air. However unlike oars in water, insects cannot lift their 
wings out of air. This leads to two subtle consequences. First, a restriction on the degree of 
asymmetry in the strokes in order to produce sufficient lift. It corresponding to an upper limit in 
the inclined angle of the stroke plane, about 60°, similar to the value observed in dragonfly flight. 
Second, a break of symmetry in the forces produced by symmetric strokes. 

PACS numbers: 87 . 19 . St ,47 . 11 . +j ,47 . 32Cc 
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Airplanes and helicopters are airborne via aerodynamic lift, not drag. However, it is not 
apriori clear that nature should design insects to fly using only lift. In rowing and swimming, 
we make use of drag to propel, so there is a reason to suspect that insects might do the 
same. 

It is well appreciated that micro-organisms (bacteria, sperm, and protozoa) use drag to 
swim in low Reynolds number flows 111 |2j. It is also known that some birds and fish use drag 
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to fly and swim in high Reynolds number flows [3, 4]. But recent research on insect flight 
has primarily focused on the , and seems to 

have overlooked the useful effects of drag. 

The separation of lift and drag, which are force components orthogonal and anti-parallel 
to the instantaneous velocity of the wing relative to the far field flow, is natural for airplane 
wings, propeller and windmill blades, and boat sails. This is because large wings and 
sails 'fly' at a relatively steady state and at small angles of attack, so lift is the dominant 
component. But an insect uses large angles of attack to generate high transient lift, i.e. 
taking advantage of dynamic stanff. At large a„g.es of attack, high lift and high drag go 
hand in hand, as expected for stalled flow and as seen in recent experiments [8 J. Therefore 
the conventional separation of lift and drag is no longer of central interest. The relevant 
force for hovering is the net vertical force which balances the gravity; thus it is convenient 
to decompose the forces into vertical and horizontal components. If a wing is restricted to 
move symmetrically along a horizontal plane, whether in circular helicopter-like or reciprocal 
motions, the drag roughly cancels in each cycle. And the net vertical force to balance an 
insect's weight only comes from lift. In contrast, asymmetric up and down strokes along 
an inclined stroke plane generate vertical forces which are oblique to the stroke plane. The 
vertical force thus has contributions from both lift and drag. 

So at least in theory, drag and lift can be of equal use to insect. But can using drag be 
as efficient as using lift at the range of Reynolds number appropriate for insects? 

Strong evidence suggesting that some insects use drag to hover goes back to an early study 
by Weis-Fogh 9] who noticed that true hover-flies (Syrphinae) and dragonflies (Odonata) em- 
ploy asymmetric strokes along an inclined stroke plane. This is in contrast to the symmetric 
back and forth strokes near a horizontal plane, seen in the majority of hovering insects, 
including fruit-flies, bees, and beetles. What brought the asymmetric strokes to attention 
was the failure of quasi-steady calculations in predicting enough vertical forces to support 
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an insect's weight 

Bay 

, 0, On the other hand, the plunging motion during the 

down-stroke in the asymmetric strokes certainly indicates that large drag can be generated 
in the upward direction. 

Since the work of Weis-Fogh, the asymmetric strokes have been treated separately from 
'normal hovering' ?| • Most studies have focused on symmetric strokes along a horizontal 
plane. However, commonly seen insect wing strokes can deviate from the horizontal plane 
and follow either parabola-like or figure-eight trajectories ^J. Here we suggest that these 
strokes can be viewed, conceptually, as pairs of asymmetric strokes, illustrated in Fig. Id and 
le. Because of this connection, we can focus on asymmetric strokes along straight inclined 
planes, and deduce qualitative results for general symmetric strokes used by insects. 

Let us first consider two cases of hovering; one employs symmetric strokes along the 
horizontal plane, a special case of normal hovering, the other asymmetric strokes alon g an 
inclined plane of 60°, suggested by the study of dragonfly free flight (Aeschna iuncea) jl0l| . 
Both cases belong to the family of wing kinematics described by 

[x(t),y(t)} = ^(l + cos27r/t)[cos/3,sin/3] (1) 
a(t) = a + Bsm(2n ft + (2) 

where [(x(t),y(t)] is the position of the center of the chord, and (3 is the inclined angle of 
the stroke plane (see Fig. la-lc). a(t) describes the chord orientation relative to the stroke 
plane. / is the frequency, Aq and B are amplitudes of translation and rotation, and is 
the phase delay between rotation and translation, ao is the mean angle of attack, thus it 
describes the asymmetry between the up and down strokes. «o = tt/2 and (3 = correspond 
to a symmetric stroke along a horizontal plane (Fig. lc). The wake of this family of wing 
motions was visualized previously using smoke trajectories jl^l]. 

The two dimensional flow around a hovering wing is governed by the Navier-Stokes equa- 
tion, which is solved here using a fourth order compact finite-difference scheme jisl]. To 
ensure sufficient resolution at the edge of the wing and efficiency in computation, elliptic 
coordinates fixed to the wing, (/i, 6), are employed and mapped to a Cartesian grid. The 
two-dimensional Navier-Stokes equation governing the vorticity in elliptic coordinates is 



«g> + (VSu • VV = ^ (3) 



3 



V • {y/Su) = 0, (4) 

where u is the velocity field, u the vorticity field, and S the scaling factor S(p, 9) = cosh 2 p — 
cos 2 9. The forces are calculated by integrating the fluid stress on the wing. Specifically, 
F = F p + F u + pWa, where F p and F^ are pressure and viscous contributions, given below, 
W is the area of ellipse, a the wing's linear acceleration, 



F p = pv J ^ (sinh p sin 9 x + cosh p cos 9y)d9, (5) 
F v = pv J lu(— cosh p sin 9 x + sinh p cos 9 y)d9. (6) 

The fictitious forces introduced by the rotating frame, i.e. the centrifugal and Coriolis forces, 
as well as the force due to rotational acceleration, integrate to zero and thus have no con- 
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tribution. The method was described in detail previously 

The instantaneous forces are nondimensionalized by 0.5pu'j ms c, where p is the density 
of air, u rms the root mean square of the translational velocity of the center of the wing, 
and c the chord, respectively. The dimensionless forces are called force coefficients, Cl and 
Co denoting the lift and drag coefficients, Cy and Ch are the vertical and horizontal force 
coefficients. Because the horizontal force cancels over a period, its absolute value is used 
when taking averages. 

The translational motion of the wing is completely specified by two dimensionless param- 
eters, the Reynolds number, Re = u max c/u = iifAoc/v, and A /c. The typical Reynolds 
number of a dragonfly is about 10 3 , and a fruit-fly is about 10 2 . But since the conclusions 
drawn below does not crucially depend on the Reynolds number in this range, it is chosen 
to be 150 for simplicity. 

In the case of symmetric stroke (Fig. lc and Fig. 2a), described by eq. (1) with a — 7r /2, 
and f3 = 0, each half-stroke generates almost equal lift in the vertical direction, and almost 
equal drag in the opposite horizontal direction. The drag in each stroke does not exactly 
cancel even in the steady state due to a symmetry breaking, as will be discussed. The 
averaged vertical and horizontal force coefficients are 1.07 and 1.61, respectively, resulting 
in a ratio of 0.66. In contrast, the asymmetric stroke with a = 60° (3 = 62.8° (Fig. 3a and 
3b) generates most of its vertical force during the down-stroke, in which the lift and drag 
coefficients are 0.45 and 2.4, respectively; they are 0.50 and 0.68 during the upstroke. The 
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vertical and horizontal force coefficients averaged over one period are 0.98 and 0.75, resulting 
in a ratio of 1.31, twice the value of the symmetric stroke. In this case 76% of the vertical 
force is contributed by the aerodynamic drag. 

Another view of the difference between the symmetric and asymmetric strokes is revealed 
by the averaged flow around the wing. The vorticity and velocity fields around the wing at 
four different times are presented in Fig. 2c and 3c. Comparing the traveling distances of 
vortex pairs in the two cases suggests a faster jet produced by the asymmetric stroke. A 
better way to quantify these jets is to plot the time averaged velocity below the wing. The 
averaged flow shows the structure of the jets, as shown in Fig. 2d and 3d. The velocity is 
plotted in physical space, which are interpretated from the computed velocity in the body 
coordinates. The asymmetric stroke generates a faster jet of a width comparable to the 
chord, and it penetrates downward for about 7 chords. In contrast, the symmetric stroke 
generates a jet whose width is comparable to the flapping amplitude, and it penetrates down 
for about 4-5 chords. Thus we tentatively conjecture that dragonflies take advantage of the 
ground effects to hover above the water at a distance of 7 chords or less, as the jet is reflected 
by the water. 

Next we investigate how the forces and power vary with the degree of asymmetry. For 
this purpose, it is convenient to fix A /c = 2.5, / = 1, Re = 150, B = ir/4, = 0, 
and vary a from it /A to ir/2, with ten a's equally spaced. The time-averaged forces and 
power are plotted against f3, the angle of the stroke plane, in Fig. 4. For a given ao, (3 is 
determined such that the net force averaged over a period is vertically up. Fig. 4 illustrates 
two interesting points. 

First, as the stroke plane tilts up, the vertical force coefficient, Cy, remains almost a 
constant up to (3 ~ 60°. The horizontal force averages to zero, but its average magnitude, 
Cjj decreases with (3. Thus, the ratio, Cy/Cn, increases by a factor of 2 as (3 increases from 
0° to 60°. Therefore by employing asymmetric strokes along an inclined plane, an insect 
not only maintain the vertical force but also reduces horizontal forces. The averaged power 
exerted by the wing to the fluid is given by P = (Fo(t)u(t)) , where Fn{t) is the force parallel 
to the translational velocity of the wing, u(t). Comparing this power with the ideal power 



based on the actuator disk theory 



17j gives a non-dimensional measure, 
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c p = \:,;;;v h, (7) 



lis category. One 
loj | : other studies 



(Fp(t)u( t)) 
(Fv(t)Y 

where the size of the actuator for a two dimensional wing is assumed to be the amplitude 
A , and Fy is the vertical force. Similar to Cy, Cp is relatively independent of (3 up to 
(3 ~ 60°. Up to 40°, there is a slight decrease in power required to balance a given weight. 
Thus using drag to hover does not require extra work. 

Second, the sharp decrease in vertical forces at (3 ~ 60° suggests that asymmetric strokes 
along nearly vertical stroke planes would not generate sufficient force to hover. A quasi- 
steady model of the forces would not predict such a cut-off ^J. The drop in vertical force 
is due to the fact that an insect cannot lift its wing out of air, so the wing must interact 
with its own wake. Such interaction reduces the net vertical force by more than half when 
(3 > 70°. Thus, one might expect that in natural flight the angle of the stroke plane also 
has an upper limit. Insects using inclined stroke planes appear to fall in t 
of the largest angles, (3 ~ 60°, is observed in dragonflies (Aeschna juncea) 
reported smaller inclined angles 

The interaction between the wing and its wake at these Reynolds numbers has another 
consequence. Namely, symmetric strokes do not necessarily generate symmetric forces, as 
would be the case at zero Reynolds number flow. At the Reynolds numbers studied here, 
the symmetric stroke along the horizontal plane (Fig. lb) produces a slightly nonzero 
horizontal force. The symmetry is broken by the initial conditions, whose effects persist in 
time. Such persistence may seem peculiar at first, but can be explained as follows. The 
first back and forth stroke generates two pairs of vortices, one at the leading edge and 
another at the trailing edge. The strengths of these vortex pairs are asymmetric due to 
the initial motion. Each newly generated vortex is paired with the previous vortex, thus 
it 'remembers' the preceeding stroke. In other words, a sequence of vortices are 'zipped' 
together. Consequently, the vorticity configuration at a given instance can be traced back 
to the first stroke, thus preserving the initial asymmetry. The nonzero averaged horizontal 
force, the pwalternating magnitude in the peaks of the lift (Fig. 2b), and the bas in the 
average flow (Fig. 2c) are manifestations of this asymmetry. Similar broken of asymmetry 
was observed in a previous theoretical study of buttefly strokes |2o| 

Real insect wing motions are complex and diverse, and the view of the flapping motion 



taken here is simplistic. In addition, the computation is two dimensional, but insects live in 
three dimensions. Nonetheless, the above analysis provides a new way of thinking about the 
wing kinematics employed by the insects. In particular, what was perceived as an exception 
in the biology literature, the asymmetric strokes along an inclined stroke plane, can be 
viewed as building blocks of more commonly seen symmetric strokes, including the well- 
known figure-of-eight strokes. Moreover, the upper limit of the stroke plane angle predicted 
from the computation, 60°, coincides with the maximal angle observed in insects 10]. It also 
turns out that at least within the family of strokes studied here, using drag is as efficient as 
using lift. We thus hope that biologists will investigate insect wing kinematics with a new 
question in mind, that is to what degree insects use drag? Our study further suggests two 
general lessons. First, to theorists, it seems to be more natural to view insects as swimmers 
in air rather than small airplanes. Thus instead of refining the unsteady lifting line theory for 
attached flow, as was appropriate for airplanes, it is more relavent to find a better theory for 
separated flow. The classical theory of separated flow underpredicts the forces j^lj]. Second, 
designing micro-scale flapping mechanisms at very low Reynolds numbers need not follow 
the traditional rule of optimizing lift, but instead could use both lift and drag. 
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FIG. 1: A family of asymmetric strokes studied in this paper. Solid lines represent chord positions 
during a stroke. The leading edge is indicated by either a black dot (down-stroke) or an open circle 
(up-stroke). The numbers next to the leading edges indicate the time sequence. The motion is 
described by eq. (1) and (2), with Aq being the end to end amplitude of the stroke, (5 the inclined 
angle of the stroke plane (dashed line), and «o the initial orientation of the wing, a) definition of 
the chord position and orientation, b) generic asymmetric strokes c) a special case of symmetric 
strokes along the horizontal plane, d) symmetric strokes along a parabolic stroke plane, and e) 
symmetric strokes along a figure-eight. The symmetric strokes in d) and e) can be decomposed into 
symmetric pairs of the asymmetric strokes described in b). The numbers next to the wings both 
indicate the time sequence during a stroke, and also identifies the segments on the right hand side 
in the original stroke. This decomposition allows one to deduce qualitative results for the general 
symmetric strokes based on the results of asymmetric strokes studied here. 
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FIG. 2: A special case of normal hovering: symmetric strokes strictly along a horizontal plane, 
a) wing kinematics given by eq. (1) and (2) with Aq/c = 2.5, «o = tt/2, B = vr/4, = 0. The 
averaged lift and drag over each stroke are shown in scale. In this case, drag almost cancels in two 
consecutive strokes, and the vertical force is contributed primarily by aerodynamic lift, b) time 
dependent vertical (Cy) and horizontal (C#) force coefficients, which are the same as lift (G£) and 
drag (Cd) coefficients, c) snapshots of vorticity (red: counterclockwise rotation, blue: clockwise 
rotation) and velocity field (red vectors) near the wing, and d) time averaged velocity field over one 
period, characterized by a downward jet. u max , the maximum translational velocity of the wing is 
the reference scale. The dashed square corresponds to the region shown in lc). 
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FIG. 3: Hovering along an inclined stroke plane using asymmetric strokes, as seen in dragonfly 

I — I 

flight 10]. a) wing kinematics given by eq. (1) and (2) with Aq/c = 2.5, ao = 7r/3, B = 7r/4, and 
j3 = 0.357T. The averaged lift and drag over each stroke are shown in scale. In this case, the vertical 
force is contributed primarily by the drag during the down-stroke, b) vertical (Cy) and horizontal 
(Ch) coefficients (solid lines), and contributions from lift (dashed line), and drag (solid line with 
circle), c)-d) see the legend of Fig 2. 
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FIG. 4: Comparison of the averaged vertical (Cy), horizontal (Cjj), and power (Cp) coefficients 
as a function of the angle of the stroke plane, 0, which characterizes the degree of asymmetry 
in the up- and down- strokes and the relative contribution of lift and drag to the vertical force. 
The vertical force remains roughly constant up to (3 = 60°, but have a sharp cut-off beyond that, 
due to the interaction between the wing and its own wake. The magnitude of the horizontal force 
decreases with (3. Cp is relatively independent of (3 up to j3 ~ 60°. Up to 40°, there is a slight 
decrease in power required to balance a given weight. Thus using drag to hover does not require 
extra work. 
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